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Here we will derive the general theory of the beam-breakup instability in recirculating linear 
accelerators with coupled beam optics and with polarized higher order dipole modes. The bunches 
do not have to be at the same RF phase during each recirculation turn. This is important for the 
description of energy recovery linacs (ERLs) where beam currents become very large and coupled 
optics are used on purpose to increase the threshold current. This theory can be used for the 
analysis of phase errors of recirculated bunches, and of errors in the optical coupling arrangement. 
It is shown how the threshold current for a given linac can be computed and a remarkable agreement 
with tracking data is demonstrated. The general formulas are then analyzed for several analytically 
solvable cases, which show: (a) Why different higher order modes (HOM) in one cavity can couple 
and cannot be considered individually, even when their frequencies are separated by much more 
than the resonance widths of the HOMs. For the Cornell ERL as an example, it is noted that 
optimum advantage of coupled optics is taken when the cavities are designed with an x-y HOM 
frequency splitting of above 50MHz. The threshold current is then far above the design current of 
this accelerator, (b) How the x-y coupling in the particle optics determines when modes can be 
considered separately, (c) That the increase of the threshold current obtainable by coupled optics 
and polarized modes diminishes roughly with the square root of the HOMs' quality factors. Therefore 
the largest advantages are achieved with cavities that are not specifically designed to minimize these 
quality factors, e.g. by means of HOM absorbers, (d) How multiple-turn recirculation interferes 
with the threshold improvements obtained with a coupled optics. Furthermore, the orbit deviations 
produced by cavity misalignments are also generalized to coupled optics. It is shown that the BBU 
instability always occurs before the orbit excursion becomes very large. 



I. INTRODUCTION 

In several applications of linear accelerators the 
charged particle beam passes through the accelerating 
structures more than once after being lead back to the 
entrance of the linac by a return loop. By this method the 
linac can either add energy to electrons several times, or 
it can recapture the energy of high energy electrons after 
they have already been used for experiments. The former 
technique is referred to as recirculating linac, the latter 
as energy recovery linac (ERL) Q]. 

ERLs have received attention in recent years since they 
have the potential to accelerator currents much larger 
than those of non-recovering linacs, and since they have 
the potential for providing emittances smaller than those 
in x-ray storage rings at similar energies and for similar 
beam currents. This is due to the fact that the emittances 
in an ERL can be as small as that of the electron source, 
if emittance increase during acceleration can be avoided. 

Several laboratories have proposed high power ERLs 
for different purposes. Designs for light production with 
different parameter sets and various applications are be- 
ing worked on by Cornell University |d, Daresbury 
0, TJNAF 0, JAERI Q, Novosibirsk 0, and KEK 
Q. TJNAF has incorporated an ERL in its design of an 
electron-ion colhder (EIC) ^ for medium energy physics, 
while BNL is working on an ERL-based electron cooler 
for the ions in the relativistic ion collider (RHIC). 
The work at TJNAF, JAERI and Novosibirsk is based 
on existing ERLs of relatively small scale. The first in- 
ternational ERL workshop with over 150 participants in 
early 2005 has also shown the large interest in ERLs that 



is prevalent in the accelerator community. 

One important limitation to the current that can be 
accelerated in ERLs or recirculating linear accelerators 
in general is the regenerative beam-breakup (BBU) in- 
stability. The size and cost of all these new accelerators 
certainly requires a very detailed understanding of this 
limitation. In [Tlj | we have described this theory for par- 
ticle motion in one degree of freedom. Here we generalize 
this theory to two degrees of freedom, i.e. to accelerators 
with polarized HOMs and x-y coupling of the particle 
optics. 

For one degree of freedom, a theory of BBU instability 
in recirculating linacs, where the energy is not recovered 
but added in each pass through the linac, was presented 
in 12]. This original theory was additionally restricted 
to scenarios where the bunches of the different turns are 
in the linac at about the same accelerating RF phase, 
such as in the so-called continuous wave (CW) operation 
where every bucket is filled. Tracking simulations [l^ 
compared well with this theory. This theory determines 
above what threshold current /th the transverse bunch 
position X displays undamped oscillations in the pres- 
ence of a higher order mode (HOM) with frequency uj\. 
If there is only one higher order mode and one recircu- 
lation turn with a recirculation time tr in the linac, the 
following formula is obtained for T12 sinujxtr < 0: 



where c is the speed of hght, e is the elementary charge, 
{R/Q)xQ\ is the impedance (in units of il) of the higher 
order mode driving the instability, Qx is its quality fac- 
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tor. In the case of one degree of freedom, Tj*2 = T12 is 
the element of the transport matrix that relates initial 
transverse momentum before and x after the recircu- 
lation loop. A corresponding formula had already been 
presented in Occasionally, additional factors are 

found when this equation is stated [THl IT^ IT^ . But in 
it has been shown that no such additional factors are 
required. 

The beam transport element T12 appears since a HOM 
produces a transverse momentum during the first pass 
of a particle. This produces a transverse position of a; = 
Ti2Px when the particle traverses the HOM for a second 
time, and this in turn excites the HOM itself by means 
of the wake field Wx = WTi2Px- 

When the mode is polarized with an angle 0, the 
kick produced during the first turn corresponds to the 
momentum {px,Py) — p{cos 6 , sin 6) . With a coupled 
optics, the resulting orbit displacement when the par- 
ticle reaches the HOM after the return loop is x — 
p{Ti2 cos 6 + Ti4 sin 6, T32 cos 6 + T34 sin 9) . This excites 
the higher order mode by the projection of this displace- 
ment onto the wake field, W -x = Wp(Ti2 cos^ 9 + (T14 -|- 
T32)sin6lcos6l-Kr34 sin^ 9). 

The HOM therefore produces a transverse kick that 
feeds back to itself, exactly as in the case with one degree 
of freedom, only that T12 needs to be replaced by 

r*2 = T12 cos^ 9 + (ri4 + T32 ) sin 61 cos 61 -I- T34 sin^ 6* . (2) 

While equation ^ is derived with one HOM, for one 
degree of freedom it is often a good approximation even 
when the cavity has several higher order modes. It was 
shown in that different HOMs can be treated individ- 
ually when their frequencies differ by more than about 
T^^. This statement does not hold for two degrees of 
freedom as will be shown in this paper. Modes cannot 
in general be treated independently, even when their fre- 
quencies are separated by much more than the width of 
the HOM's resonance. 

An optics configuration that makes T12 close to zero in 
order to make the threshold current very large has been 
proposed flol, |2^ and tested for coupled beam transport 
and polarized HOMs. This is a good technique when 
there is one dominant HOM. When there are several 
modes this approach does not apply directly, even if these 
modes are separated by more than the width of their res- 
onance. 

The paper is arranged as following: first, a disper- 
sion relation for the current Io{lo) is derived including 
coupling in a one turn recirculating linac with one cav- 
ity having multiple polarized HOMs. The smallest real 
value of Iq that can be obtained with real lo determines 
the threshold current. Analytical solutions are given for 
the case of two polarized HOMs in one cavity. It is ex- 
plained how the dispersion relation for this simple case 
can be solved efficiently on a computer, and comparisons 
to analytical approximations are presented. Approxima- 
tions are then given for N polarized modes in one cavity. 
Subsequently a dispersion relation for multiple cavities 



and multiple recirculation loops is derived that can only 
be solved numerically with similarly efficient techniques. 
Finally, misalignments of cavities are considered to in- 
vestigate when these misalignments lead to a very large 
static displacement of the beam orbit. 

II. iV POLARIZED MODES IN ONE CAVITY 

For simplicity we are here investigating one cavity with 
N higher-order dipole modes (HOMs), each having a po- 
larization angle 9\ to the horizontal. Note that slightly 
polarized cavities often have at least two HOMs with 
similar characteristics, and looking at a single polarized 
HOM can therefore be misleading. 

The unit vector in the direction of the polarization is 
ex = cos9\ex +sn\9\ey. The effective transverse voltage 
in each HOM is V\, so that a particle traversing this 
HOM obtains a transverse-momentum change of ^e\V\. 

When V is the vector of all these N voltages, then the 
momentum change is 

Ap= -EV , E = (f?i...ejv) . (3) 

c ^ ' 

When the particle returns to the cavity after the return 
time tr, the particle's position has changed by 

Ax{t') = TAp{t' - U) ^ T^=(^tZtH) ■ ("^^ 

This position will increase the voltage in Vx{t) by the 
projection of the mode's wake field W\{t — t')e\ onto 
Ax{t') times the charge I{t')dt' that excites the field, 
i.e. I{t')Wx{t - t')ex ■ Ax{t')dt'. Integrating over all 
contributions to the HOM potentials leads to 

V(t) = - f I(t')W(t ~ t')T>V(t' - tr)dt' (5) 

C J-00 

with 

W = d\a.g{Wx) , D = T E (6) 
and therefore 

= T12 cos^ 9^ + sin2 9^ + ^ sin 29^ , 

Dfj^i, = T12 cos 9^ cos 9 1, + r34 sin 9^ sin 9i, 

+ r32 sin 6*^ cos 9u + T14 cos 9^ sin 9,^ . (7) 

Here I{t') is the current of the bunches that have al- 
ready traveled for one turn; in the approximation of short 
bunches it is given by 

00 

I{t) = Iotb 5{t-tr-mtb). (8) 

This transforms the integral equation into 

00 

V{t) = Iq-U ^ W(t -U- mtb)'DV{mtb) (9) 
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where Wx{t) = for i < 0. 

The Laplace transform of V{t) can be written as 

Vit + Stb) = — / VsiLj')e''^''du;' . (10) 

And with the foUowing definition 



27r 



(11) 



n— — oo 



one obtains 

y,^H=4 ^ y([n + J]4)e^'^"*^ . (12) 



n— — oo 



Since V/'((jj) is periodic with 27r/th, it has a Fourier 
series, and its Fourier coefhcients are V^([n + (JJtb), which 
shows that V/'(w) does not vanish. The transverse mo- 
tion is stable when V{ijj) is zero for all uj with positive 
imaginary part. If the current is increased the motion 

can become unstable at which point V{uj) is non-zero for 
at least one w with positive imaginary part. At threshold 
it is therefore non-zero for a real value of uj. 

As in we will use tr = {rir — S)tt, to allow for all 
recirculating phases, e.g. 5 = ^ for an ERL. The integral 
equation now leads to a relation for these coefhcients. 



(13) 

W([n -nr-m + 5]tb)T>V{mtb)e"^''*'' 



m.n— — OO 



= -e^^^-'-'Wf (w)Dt/cp(w) . 

This formulation shows that V^{uj) is an eigenvector of 
the matrix on the right-hand side, and the corresponding 
eigenvalue is I/Iq. Its solution is therefore very similar to 
the matrix theory for BBU computation without coupling 
in p^. The threshold current can thus be determined 
by finding the largest real eigenvalue of this matrix for 
any real w. Due to the symmetry properties W|"(a; -|- 
1^) = Wf (w) and Wf{-Lj) = Wf*(w), it is sufficient 
to investigate uj e [0, 7r/ib] to find the BBU threshold 
current, 

/^-i = 1 max {A|e*'^"'-*''Wf (tj)Df = AV,AeM}. 



A. Two polarized HOMs in one cavity 



(14) 



For a large number of HOMs this equation should be 
solved numerically, but for two HOMs a analytical so- 
lution is simple. The characteristic polynomial for the 
eigenvalue I/Iq becomes 



/ ^ - Diiwi -D12W1 
-D21W2 - D22W2 



, (15) 



with / = ^ike^^'-Jo and wx = ^e*"**" (Wf )aa. Solving 
this quadratic equation leads to 



c 

etb 



e "in = 7: 



(16) 



-DllWl - D22W2 y2 
±\j( Y + WiW2-L>i2D21 



To reduce the threshold current, the left hand side of this 
formula should be small. The matrix D is determined by 
the mode polarization and by the linear particle optics. It 
has been suggested 0, |3 to use these parameters to in- 
crease the threshold current. These suggestions amount 
to reducing the right hand side of Eq. H16|l by making 
Dii and D22 small. This is always a valid strategy if 
there is only one HOM, i.e. W2 = 0. Then 



Dii = T12 cos^ 6*1 + T34 sin^ Oi 



32 



sin 26*1 (17) 



and a horizontally {9i = 0) or a vertically {di — ^) 
polarized mode together with a beam transport that fully 
couples the vertical to the horizontal motion and vice 
versa, i.e. T12 — T34 — 0, would always lead to Dn — 
so that there would not be any threshold current. A 
formula for this case has been derived in . 

For the case of two or more HOMs this method is no 
longer as effective, even when the modes have very dif- 
ferent frequencies. In fact it has been suggested that in 
cases with several modes, each mode could be considered 
separately when each cavity mode has a resonance width 
which is significantly smaller than the frequency sepa- 
ration between modes This is however not correct 
and it will be seen shortly that the described method 
that would seem to increase the threshold current when 
all modes are considered separately is not as effective 
even when modes frequencies differ by much more than 
the resonance width. This is especially important when 
there are two HOMs of similar properties, for example in 
nearly cylindrically symmetric cavities. To see this effect, 
we distinguish three cases. 



1. Circular symmetry 

For circular symmetric cavities there are two equivalent 
modes with perpendicular polarization, wi = W2, ei — 
ex, £2 = Sy leading to D = T, 



/n = 



ctb T±wi 



(18) 



T±e ± = ±y( y + TiiT32, 

with T±eIR+ . 

where the threshold current is the smaller of the two val- 
ues obtained with the equation for -I- and for — . For 
T12 = T34 = this case has been considered in [23|. If 
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there is no coupling, T+ = \Ti2\ and T_ = IT34I. This 
result is equivalent to what has been found in ITl| for 
motion in one degree of freedom. 

As in now use long range wake fields of the 

form 



Wx{t) = ^ ^e"^ sinKr) , (19) 

where Q\ is the quality factor, wx is the frequency, and 
{R/Q)x is the impedance in units of Vl for the linac defi- 
nition (2 times the circuit definition). As shown in [Tll |. 
|z«a(w)| is especially large when lo is close to 

When motion in only one dimension is considered as 
in one obtains 



c e 



etb T12WX 



(20) 



For simplicity we again use K.\ = h-j^ (^^^ and ex = 



For UrCx 1 this leads to the approximation 
' 1 forri2sin(wA) <0, 



Ith 



Kx T12 sin(a)Atr) 
2 



+ ( 



mod(i£AW[) ^2 



)2 else. 



(21) 



I k;a|Ti2| ■ 

For ^ 1 but nr-ex 3> 1 the approximation derived in 
[Til is 



Ith 



2eA 



(22) 



ICx\Ti2. 

Exactly the same approximation and derivation there- 
fore leads leads from Eq. 118|l to 



Ith± 



(23) 



For ea ^ 1 but nr€x 3> 1 the approximation derived in 
ini is 



th± 



2ex 
ICxT± 



(24) 



The threshold current is given by Ith = min±(/tft±). 

To clarify the here considered case, we use typical 
parameters for the two HOMs: Qi = Q2 ^ 10'*, 
(i) = (i) = ^'^'^^^ uji = LU2 = 2tt ■ 2.2GHz, 
tb = l/1.3GHz, rir — (5 = 5.5. For a decoupled optics with 
T12 = -W-^-^, Tii = T32^0 we obtain a threshold 
current of Ith = 46.40mA that agrees to all specified dig- 
its when computed by particle tracking and by the ap- 
proximation in Eq. 124|l . For a very much coupled beam 
transport (abbreviation x = ^IO^^^^tt^ used below) 
the following threshold current is obtained: 



1 X 1 3a; 

10 1 

-1 -2x 1 Ax 
0-101 



ith 



20.28mA 



(25) 

Again the threshold current computed by particle track- 
ing and by Eq. (|23|l agree to all specified digits. 



2. Small coupling 

To cases with \DiiWi — -022^^2 P ^ |wiU'2£'i2-D2i| we 
refer to it as small coupling. This denomination is moti- 
vated by the fact that when one mode is polarized in x 
and one in y direction, this case occurs when the optical 
coupling is small, as can be seen in Eq. Q. Equation 
(|16|l simplifies to 



ctb Dxxwx 



(26) 



with A being either 1 or 2. 

The approximation that leads from Eq. (|20|l to Eq. 121|) 

can again be used and it leads to 



2eA 



th 



tCx DxX SV!l{u>xtr 



ior DxxsH^xU) <0 , (27) 



and similarly an approximation that corresponds to 
Eq. H22I) is valid. This formula has been used to argue 
that an optics with very small |-Daa| could be built with 
extremely large threshold current. But this equation does 
not apply when the |-Daa| are too small, since then the 
following case has to be considered. 



3. Strong coupling 

To the case with \DiiWi — D22W2P ^ IW1W2D12-D21I 
we refer to as strong coupling. This case is especially 
relevant when the mentioned coupling techniques have 
been used to make Dn and D22 very small. One obtains 



-iuj2tj. 



{—f— 

eth W1W2D12D21 



(28) 



At the threshold current only one of the functions in the 
denominator will be very large and we call this w^. The 
other mode will be indexed by v. 

In jlH a first order approximation in Aujtb = {ui—ujf^)tb 
and is used to obtain Eqs. 1)21(1 and 1)22(1 . Here we again 
expand to first order in these quantities, which is simple 
since is linear in them so that can be evaluated 
at Acj = and = 0, leading to 



c e 



-i(u2tr+^) 



etb Twi 



(29) 



With Te'^ = ^w^{uj^)Di2D2i, T € IR+ . Note that the 
exponent contains 2tr instead of tr . The approximations 
that correspond directly to Eqs. H23() and (|24() can again 
be applied. 

In hj wx is evaluated for the long-range wakefield, 
leading to 



R 

e^^+iS-i)tbsi-a{LJxStb) 



(30) 



cos 



'^'^^*''sin(a;A[^- 1]4) 
'.{uj+tb) - cos{ujxtb) 
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Note that here wx is defined without a phase factor 
^-luStb compared to to simphfy the notation. For 
an ERL, where 6 = 1/2, one has 



Wx = 



\Q J X 2c cos UJ+tb — cos LUxtb 



OJ 



.2 co s(a;+|)sin(a;A|) ^^^^ 



with w+tb cijtft + ieA- Then lu — lu^ + Auj with the smaU 
quantity Auiti,. We assume that also = is small, 
which is usually the case whenever BBU is relevant. A 
first order expansion in these small quantities leads to 



2e^ 



'"^(Awtb + iSfj,) cos(tj^tb) - cos(wytf,) 



K.U.IC1JD12D21 



cos(tj^^)sin(tJi.^) 



(32) 



If /q is the threshold current Itin the right hand side 
has to be a real number, requiring Auitb sin(ct'2tr) = 
ei cos(ti;2tr)- This leads to 



2e„ 



COSLOntb — COS (jj^tb 



thfi 



JC^JC^ D12D21 cos(w^,^) sin(wi.^) sin(a;^2tr 



(33) 

whenever this term is positive. Whenever it is negative, 
the following approximation follows from [Tll |. 



7^2 



COS UJn,tb — COS LO^tb 



D12D21 cos(w^f ) sin(w^f ) 



(34) 



1^2 ^ ^ mod (^^2tr,7r) ^^ 



2nr 



For ^ 1 but n^^p 3> 1 one obtains 



7^2 



2e„ 



COSUJ^tb — COS UJ^tb 



D12D21 cos{ujf^i^) sin(a;,,^) 



(35) 



These formulas are to be evaluated for /z = 1, v = 2 
and for fi = 2, v = 1, and the smaller of the two resulting 
currents is the threshold current. 



Ith = min {Ithfj.} 



(36) 



An interesting observation is that for two modes with 
similar Qx, , and cox, the ratio of the threshold 

current with and without coupled optics can be found by 
comparing Eqs. H33|l and (|27|l and it is proportional to 
y^l/ex oc \/Qx- For cavities that are optimized for large 
currents by means of sophisticated HOM damping, the 
advantage of a coupled optics therefore decreases. This 
effect is independent of the length of the return loop and 
is already relevant for a single cavity. 

Figure ^ refers to a linac with one cavity and tr = 
555.5tfc that has two HOMs, one with = 2.2GHz po- 
larized in X and one with 102 — 2.3GHz polarized in y 
direction, = lOOil, and Q^, — Qy is varied. The 

bottom data refer to a decoupled optics and the top data 
to a fully coupled optics with T12 = and T34 = 0. A 



logio(/«,(A)) 
2 






10 



12 logio(Q) 



FIG. 1: The threshold with two modes in one cavity scales 
with for a decoupled and with Q~2 for a fully coupled 
optics. The advantage of coupling thus diminishes for low Q. 



double-logarithmic plot is shown, which makes it appar- 
ent that the threshold current for a decoupled optics de- 
creases with Q^^, as indicated by the line with slope —1. 
A closer look shows that the slope is only accurately — 1 
for relatively small and relatively large Q where either 
approximation [2] or [53 hold. A totally coupled optics 
with T12 = and T34 = leads to a larger threshold cur- 
rent than without coupling, but when Q of the modes is 
reduced, the threshold current increases only with 
as indicated by a line with slope —\. The advantage of 
coupling decreases proportionally to ^/Q. 

Figure |21 similarly shows the advantage of polarizing 
the higher order modes in the ERL that is proposed to 
upgrade the CESR ring at Cornell University [23 ■ When 
the HOMs are polarized in x and y direction and the 
optics is completely coupled by T12 = and T34 = 0, 
the threshold current is larger than without coupling, 
but again the advantage is smaller when the HOMs are 
damped more strongly by HOM absorbers. In fact, for 
the HOMs that are computed for 7-cell cavities 0| , the 
increase of the threshold current due to coupling is only 
a factor of 2.5 as can be seen in Tab when the HOM 
frequencies in x and y are separated by lOMHz. However, 
the figure shows that when many cavities are present, as 
in the ERL where there are 320, the scaling is not as 
simple as in the case of a single cavity and the advantage 
of coupling does not decrease as strongly with decreasing 
Q. 

For Tab.[T|the 320 cavities of the ERL upgrade of CESR 
had nominal HOM frequencies of fx — 1. 87394GHz with 
horizontal polarization. The mode with vertical polariza- 
tion is fy — fx — Afxy The cavities have HOM frequen- 
cies that have a Gaussian distribution around these val- 
ues with rms width ajip- We used 500 different random 
distributions of the frequencies and display the average 
threshold current Ith as well as the rms ct/ of the 500 re- 
sulting thresholds. All modes have = 109.60r2 and 
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FIG. 2: The threshold increase due to couphng of the x-ray 
ERL upgrade of CESR scales roughly with y/Q. 
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FIG. 3: Two eigenvalues trace out loops for Jo in the complex 
plain while uj is varied. 



TABLE I: Threshold currents for the two most significant 
HOMs of the Cornell ERL. 
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B. Comments about numerical solutions 

When the threshold current for two polarized modes 
at uji, and should be found by solving Eqs. (|26|l and 
(|28|l numerically, the eigenvalues are plotted in the com- 
plex plain, and the intersections with the real axis are 
sought that lead to the smallest current. An example 
of the two eigenvalues plotted in the complex plain is 
shown in Fig.|31 The eigenvalues are largest in the vicin- 
ity of UJ — \mod±{Lj\,j^)\, A g {1,2}, since there ei- 
ther or become very large. The subscript on 
the mod function indicates that mod+(a:, 1) e [0, 1] and 
mod_(a;,l) e [-^,5]. 

Furthermore, the eigenvalues trace out loops around 
the origin of the complex plain about once per 7^ varia- 
tion of Lu, due to the exponential factor in Eq. (|26|l . We 
therefore vary ui only in a interval around each HOM 
frequency. This speeds up tfie search for eigenvalues by 
a factor proportional to n^, which can be very large. 

A simple approach would be to plot all eigenvalues 
in the complex plain and to select the smallest eigen- 
value that is reasonably close to the complex plain. Large 



FIG. 4: Variation of the relative accuracy AL for scanning 
io — Lu^ G [~T^i T^] t'y ^ points. By far the best accuracy is 
achieved with the elliptical extrapolation of 4 points around 
the real axis (black dots). For an approximation by a line 
between two points (red, largest at large k), a second order 
polynomial y{x) fitted to three points (green, second largest 
at large k), and a third order polynomial y{x) fitted to four 
pints (blue, third largest at large k) all lead to much worth 
accuracies and/or computation times. 



factors in speed can be gained when the loops that are 
traced out by each eigenvector can be interpolated and 
their intersection with the real axis can be found, since 
the loops have to be scanned much less densely. We have 
found that using values of Io{uj), where two are above and 
two below the real axis, and fitting an upright ellipse to 
these values is a very good parametrization. The accu- 
racy achieve for distributing k particles in the described 
region around each HOM frequency lead to the accuracy 
documented in Fig. 01 
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FIG. 5: Threshold current for one horizontal HOM with fx ~ 
2.2GHz as a function of the frequency fy of a vertical HOM. 
Black curve: dispersion relation, Red dots: tracking, Colored 
curves: the different applicable approximations derived above. 



C. Comparison of results for polarized modes and 
coupling 

To demonstrate the excellent agreement between these 
numerical solutions and tracking, we depict Fig. [S] where 
one HOM with horizontal polarization has been fixed at 
fx = 2.2GHz, and another has been varied for fy e 
[1.2, 3.2]GHz. The optics was completely coupled with 
Ti2 = T34 = and Tgi = -T^ = ^-^WTc- ^^^^^^^^^ 
things can be observed: (1) since tracking is relatively 
time consuming, only relatively few frequencies for the 
second HOM have been evaluated, but all of them lie ex- 
actly on the curve that follows form the dispersion rela- 
tion. (2) The threshold current varies strongly when the 
second HOM is varied, but the agreement with tracking 
shows that this is not numerical noise but a consequence 
of the coupling between the two polarized modes whose 
frequencies are very far apart. (3) There are frequen- 
cies where the threshold current is relatively small, these 
are frequencies where cos{uJitb) ~ cos{u}2tb) in agree- 
ment with Eq. H32I) . The displayed minima appear at 
(3 • 1.3 - 2.2)GHz, 2.2GHz, and (4 • 1.3 - 2.2)GHz. But 
the regions with reduced threshold are relatively wide; in 
particularly they are much wider than the width of HOM 
resonances. 

Depicted in color are the values obtained by approxi- 
mations derived above. It is apparent that the approx- 
imations are not always very good, especially that they 
lead to values that are too large. A magnification in FigEl 
shows that the reason for that is that the approximate 
formulas lead to parabolic shapes that have the correct 
minimum value, but not the correct width. This is im- 
portant since it shows that the formulas can be used to 
find the correct minimum value as a conservative esti- 
mate. Furthermore, the magnifications again show the 
good agreement with tracking results. 
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FIG. 6: Magnification of regions in Fig.|3 



Here an important note is in place. Once place where a 
strong dip occurs is when uji ~ lu2. And this dip is much 
wider than the width of the HOM resonance of u!\/QX, 
showing that the two modes clearly do not decouple when 
they are separated by more than their width. For nom- 
inally circular symmetric cavities, HOMs are not degen- 
erate due to construction errors and each mode splits 
into two mode with typically a few MeV distance. But 
the dip is much wider than that, showing that an appro- 
priate advantage of BBU suppression by damping can 
in general only be realized when polarized cavities are 
designed, i.e. cavities where the horizontal and vertical 
dimensions are designed to be slightly different, leading 
to HOM frequencies that differ by several lOMeV in the 
two planes. 

The question arises how far the HOM frequencies have 
to be apart. In Fig. [7| we show for the Cornell ERL 
how Ith changes with Af^y In order to avoid averaging 
over many different frequency distributions we have here 
chosen the same two HOM frequencies for each of the 320 
cavities. The data indicates that a mode separation of 
60MHz is sufficient. This has been used to compute the 
very large recirculating BBU threshold current of close to 
2.5A for this accelerator. It should be noted that only the 
two dominant HOMs were considered in this simulation 
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FIG. 7: Dependence of Ith on Afxy for the Cornell ERL. 



and a detailed study is in order when such high currents 
are sought. The Cornell ERL is designed for 100mA and 
this calculation indicates that polarized cavities with a 
coupled optics provide a very comfortable safety margin 
with respect to this instability. 



D. Approximation for N polarized HOMS in one 
cavity 

As a first approximation one can assume that one com- 
ponent of the eigenvectors in Eq. 1141 will be very large. 
This would lead to a decoupling of HOMs so that each 
HOM could be treated separately, as for a single degree of 
freedom. One could therefore derive the threshold with 
the smallest current for each individual HOM, and this 
would approximate the threshold for the complete accel- 
erator. When all these thresholds are very large, one has 
to investigate the next approximation, where the eigen- 
vector has two dominant components. In this case the 
eigenvalues are determined from a 2 x 2 matrix corre- 
sponding to Eq. IjlSI) . The threshold current has to be 
computed for each pair of HOMs. The smallest current 
obtained for one pair of polarized modes then approxi- 
mates the threshold current of the full accelerator. 



III. POLARIZED HOMS IN MANY CAVITIES 
AND FOR MULTIPLE TURNS 



Recirculating linacs with many cavities and several re- 
circulation loops have been considered early on [l^ |23| . 
In |ll| a description for arbitrary recirculation times has 
been presented. Here we want to extend this description 
to include orbit coupling and polarized modes. As far 
as possible we retain the notation of these earlier papers. 
The N higher order modes, which can be associated with 
different cavities, are numbered by an index i. The Np 



passes through the linac are numbered by an index /. 
The horizontal and vertical phase space coordinates that 
the beam has at time t in the HOM i during turn / is 
denoted z({t) = {xj {t) , piS) , vf (t) , piS)V ■ The 4 x 4 
transport matrix that transports the phase space vector 
at HOM j during turn J to zf is denoted Tij"^ and the 
time it takes to transport a particle from the beginning 
of the first turn to HOM i during turn / is denoted tj. 
The beam is propagated from after HOM i — 1 to after 
HOM i by 



All 



,-zl,it-[tl-ti_,]) + -V,it) , (37) 



with Vi — Vi{0, cos 9i, 0, sin 9i)^ . This equation can be it- 
erated to obtain the phase space coordinates as a function 
of the HOM strength that creates the orbit oscillations. 
With the matrix 



(38) 



one obtains 
Nuii - 1) 



/ Nij{i~l) 

j=i j=i 



N, if / 7^ J; 
if/ = J. 



(40) 



The strength Vi{t) of the HOM i with polarization di- 
rection Cj is created by all particles that have traveled 
through that HOM via the integral 



W,{t~t')ll{t')ejxl{t')dt' , (41) 



1=1 



where // {t) is the current at time t that the fraction of the 
beam has which passes the HOM i on turn /. Note that 
W{t -t') = for f > t. Combining this with Eq. ^ 
leads to the following integral-difference equation: 



V,{t) = / 5]M/,(i-0//(t') 



(42) 



1=1 

I Nu(i-l) 



^ ^ Df^V,(t'-\ti-t^])dt' 



j=i j=i 



(43) 



This equation is identical to that obtained for one de- 
gree of freedom, only that (Tf^)i2 is replaced by Dfj. 
The following treatment for obtaining the threshold cur- 
rent is therefore identical to that in [ij. For completion 
it is here presented in simplified form, and recommenda- 
tions for numerical solutions are given. 

Now the approximation of short bunches is used. The 
current is given at time t by pulses that are equally spaced 
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with the distance tt, 

oo 
m— — oo 

This reduces the integral to a sum, 



(44) 



oo "p 



V^{t) - -/otb ^ ^W^,(i-if-mi,) (45) 

m — — OO /—I 

^ E E DljV,{mh + tj) . 



This leads to 



oo oo 



^Mf/t.H = -^oi^ E E E^»("*^''+*^-*') (46) 

n— — oo m— oo I—l 

X E E ^^'VS([n-m]4 + t/)e— . 

Np I Nij(i-l) 

= ^E^f-*iE E ^^^//.M- 

7=1 J=l j = l 



If a vector F is introduced that has the coefficients 
V^^j , this equation can be written in matrix form, 



— V = M{uj)V 

^0 



(47) 



with the matrix coefficients 



Np 

|i, ifj>«; 

I 0, otherwise. 



(48) 



Note that W,^,._,.]/t. = ^^s.n.Lf 



iijTop{- 



-)tb 



where 

Top(a;) is the smallest integer that is equal to or larger 
than X and Si{I,L) = mod(if — tf',ti,). With Kronecker 

Sik this determines the matrices W and U to be 



-ti,Wi{5{I, L))e 
c 



(49) 
(50) 



For each frequency uj, 1^^ is an eigenvalue of M(ti;). 
Since the eigenvalues are in general complex, but Iq has 
to be real, the threshold current is determined by the 
largest real eigenvalue of M((jj). The matrix has the 
properties 



M(w- 



2^. 
tb' 



M(w) , M(-w*) =M*(w) 



(51) 



and it is therefore again sufficient to investigate uj G 
[Ojir/tb] to find the threshold current. 

Note that V^'' never appears since the last kick on 
the last turn does not feed back to any HOM, so that 
the dimension of M can be reduced by one to N ■ Np — 
1. Furthermore the dimension can be reduced when two 
fractional parts S( and (5/ are equal since then V/ and 
V/ are identical. 



A. Multi turn operation and cavity misalignments 

Since the formalism presented here that includes po- 
larized modes and coupled optics is identical to the for- 
malism in one degree of freedom, only that (Tf^)i2 has 
been replaced by D(j , all conclusions about multi turn 
recirculating and multi turn ERLs hold. For example, 
in it was concluded that for one HOM and for Np 
passes through the linac, the threshold current should 
roughly scale as Np{2Np — 1). The origin of this conclu- 
sion results from the double sum X]j=i- Since the 
same summation appears, this conclusion holds also for 
polarized modes with coupling. 

For misaligned cavities, HOMs are excited even when 
the current is smaller than the threshold current. This 
can lead to large beam excursions, and in it was 
analyzed for what currents these excursions become ex- 
tremely large. It was found that the BBU threshold cur- 
rent is always smaller than the current for which these 
orbit excursions would get very large. Since the here pre- 
sented formalism with coupling and polarized modes has 
the same formal structure, this conclusion again holds. 



B. Comments about numerical solutions 

As pointed out above where numerical solutions for 
one return loop and two HOMs were found, it is very 
essential to systematically search for real values for the 
eigenvalues of M. Each eigenvalue traces out curves in 
the complex plain when uj is varied in the region [0, 7r/tb), 
however eigenvalue finders usually do not return eigenval- 
ues in any particular order, so that these curves cannot 
be observed easily. If they could be observed, then el- 
lipses could be fitted to these curves and the intersection 
of the curve with the real axis could very efficiently be 
found for each eigenvector. 

We therefore recommend a sorting algorithm that sorts 
eigenvalues rather robustly: (1) Normalize each eigenvec- 
tor. (2) Sort these vectors according to their largest com- 
ponent, i.e. the vector which has its largest component 
in position 1 is the first vector, if there are more than one 
of this kind, the one with the largest coefficient can be 
chosen as first vector, etc. (3) Associate the eigenvalues 
in the order of these eigenvectors. Small changes of uj do 
not change the relative size of the eigenvector elements 
much. The intersection of the curve Xi{uj) with the real 
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axis can now be found for each eigenvector. This proce- choosing the largest eigenvalue that is reasonably close 
dure leads to an enormous speed advantage over simply to the real axis, 
scanning all eigenvalues for a mesh of w e [0,7r/t(,], and 
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